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ABSTRACT 

This  paper  presents  constructive,  network-based  proofs  for  the  Wagner 
Whitin  property  and  a  generalized  Nested  Schedule  property  for  optimal  pro 
duction  schedules  in  acyclic  hierarchical  multiperiod  production  systems. 
Algorithms  for  obtaining  schedules  with  these  properties  and  described  and 
illustrated  with  examples. 
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INTRODUCTION 


This  paper  addresses  the  problem  of  lot-sizing  in  hierarchical  pro¬ 
duction  systems.  Specifically,  we  use  a  constructive,  network  based  method 
to  establish  the  Wagner-Whiten  (WW)  property  and  a  generalization  of  the 
nested  schedules  (NS)  property  for  optimal  production  schedules  in  unca¬ 
pacitated  acyclic  multistage  systems.  In  establishing  these  properties, 
we  use  a  generalized  network  flow  model  similar  to  that  of  Steinberg  and 
Napier  [4]. 

In  the  problem  we  address  here,  there  are  N  production  stages  requiring 
production  decisions  for  each  of  T  periods  in  the  planning  horizon.  For 
any  stage  i,  P(i)  in  the  index  set  for  predecessor  stages  and  S(i)  is  the 
index  set  for  successor  stages.  If  S(i)  =  then  stage  i  corresponds  to 
a  final  product,  and  if  P(i)  -  $,  then  stage  i  corresponds  to  the  initial 
processing  of  raw  material.  Define  F  -  (i|P(i)  =  $},  the  set  of  initial 
(or  first)  production  stages  and  L  ■  { i | S (i)  *  $},  the  set  of  final  (or 
last)  stages.  The  production  of  one  unit  at  stage  j,  j  i  F,  requires 
units  from  stage  1,  1  eP(j)  where  is  a  positive  integer.  For  stages 
1  e  L,  the  demand  in  period  t  is  D^t. 

The  decisions  required  are  denoted  by  X^t ,  representing  the  number  of 
units  to  produce  at  stage  i  in  period  t.  The  end  of  period  inventories 
are  denoted  by  1^.  The  costs  associated  with  production  and  inventory  are 
assumed  to  be  concave  function  and  are  denoted  by  ^it(Xit)  and  Hit(Iifc) 
respectively.  Note  that  C^fc(X^t)  can  incorporate  a  fixed,  or  setup  cost. 

We  adopt  several  standard  assumptions,  namely: 

(1)  initial  inventories  and  production  lead  times  are  zero  at 


each  stage; 


(2)  transfer  tines  between  stages  are  zero; 

(3)  backlogging  or  lost  sales  are  not  permitted;  and 

(4)  there  are  no  capacity  limits  on  production  or  inventory. 

Also,  for  convenience,  we  assume  that  if  j  e  S(i),  then  j  >  i.  Let  P  de¬ 
note  any  problem  having  the  structure  just  described. 

THE  NETWORK  MODEL 

Figure  1  illustrates  the  logical  structure  of  a  five  stage  acyclic 
problem  and  Figure  2  illustrates  the  corresponding  network  model  for  a 
four  period  horizon.  In  the  network  there  are  two  nodes  for  each  stage- 
period  combination.  The  first,  referred  to  as  a  collector  node,  corresponds 
to  obtaining  units  from  the  predecessor  stages,  while  the  second,  referred 
to  as  a  distributor  node,  corresponds  to  disposition  of  current  production:  - 
and  the  Inventory  from  the  previous  period.  Production  arcs  join  collector 
nodes  to  distributor  nodes,  and  their  flows  correspond  to  X^.  Inventory  - 

arcs  join  distributor  nodes,  and  their  flows  correspond  to  Iifc.  Transfer _ 

arcs  join  distributor  nodes  to  collector  nodes,  and  their  flows  correspond 

to  Kj^it,  j  e  P(i). 

The  nodes  in  the  flow  network  are  numbered  in  a  particular  order,  which 
is  used  later  in  the  proofs.  At  each  stage,  all  the  collector  nodes  have 
consecutive  indices,  and  all  the  distributor  nodes  have  consecutive  indices 
with  smaller  indices  corresponding  to  earlier  time  periods.  Moreover,  a 
production  -arc  starting  at  collector  node  i  will  terminate  at  distributor 
node  i  +  T.  The  first  collector  node  for  stage  n  will  have  index  2(n-l)T  +  1 
and  the  last  distributor  node  will  have  index  2nT. 
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•  Figure  1.  A  4-Stage  Acyclic  Structure 
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Figure  2.  Network  model  for  the  4-stage  acyclic  system 
of  Figure  1,  with  a  planning  horizon  of  4  periods. 


For  arc  (m,  n)  in  the  network,  the  notation  is: 

X(m,  n)  ■  flow  on  the  arc 
K(m,  n)  ■  arc  multiplier 
C(m  aj(X(m,  n)  *  cost  of  arc  flow. 

In  addition,  for  node  m  in  the  network, 

i(m)  "  the  stage  corresponding  to  node  m 
t(m)  *  the  period  corresponding  to  node  m. 

A  flow  in  the  network  is  feasible  provided: 

(1)  at  each  distributor  node,  total  incoming  flow  is  equal  to 
total  outgoing  flow;  and 

(2)  at  each  collector  node,  the  incoming  flow  on  an  arc  is 
equal  to  the  arc  multiplier  (K^)  multiplied  by  the 
outgoing  flow. 

There  is  an  obvious  correspondence  between  feasible  flows  in  the  network 
and  feasible  production  schedules.  For  any  problem,  P,  let  N(P)  be  the 
corresponding  network. 

Steinberg  and  Napier  [4]  first  utilized  a  network  model  like  this  as 
the  basis  for  a  HIP  formulation  of  the  problem.  Subsequently,  McClain  et  al., 

[3]  demonstrated  a  more  parsimonious  mathematical  formulation,  but  did  not 
perform  any  computational  assessment.  Our  objective  is  to  use  the  network 
N(P)  to>  establish  two  conditions  that  are  satisfied  by  some  optimal  solution  - 
to  P. 
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THE  WAGNER-WHITIN  PROPERTY 


The  Wagner-Whitin  property  was  first  established  for  single-stage 
problems  [7J,  and  later  extended  to  acyclic  production  systems  by  Zangwill 
[8]  and  Veinott  {5} .  In  essence,  it  says  there  is  at  least  one  optimal 
solution  for  which  production  and  incoming  inventory  are  not  both  positive 
in  any  stage-period.  Yet  another  proof  for  the  property  is  offered  here 
because:  (1)  our  proof  is  constructive  and  network-based,  thus  more  in¬ 
tuitive  than  the  analytic  proofs  of  [5]  and  [8];  and  (2)  the  flow  adjustment 
operation  used  in  the  proof  is  also  required  in  our  proof  of  a  generalized 
nested  schedules  property. 

Theorem  1.  There  exists  an  optimal  solution  to  any  problem  satisfying 

XitIi,t-l  0  for  *  =  1 . N 

and  t  -  1,...,  T 

Theorem  1  can  be  restated  in  terms  of  the  network  M(P)  as  follows: 

Theorem  IN-:  -  There  exists  an  optimal  flow- in  N<P)  such  that  at  most  one  —  - 
arc  entering  each  distributor  node  has  a  non-zero  flow. 

We  shall  prove  the  network  version  of  the  theorem.  Our  method  of  proof  is 
similar  to  that  used  by  Wagner  (7],  i.e. ,  we  assume  that  an  optimal  flow  is 
given  whichviolates-the  theorem  and  show  how  to  modify  it  to  obtain  a  flow  . 
satisfying  the  theorem  without  increasing  the  cost. 

To  start  with,  we  have  a  given  feasible  flow  that  does  not  satisfy  the 
WW  property.  We  consider  the  distributor  nodes  in  the  order  of  Increasing 
index,  and  stop  with  the  first  one  violating  the  WW  property.  For  this  node, 
let  the  incoming  flow  on  the  production  arc  be  U  units  and  that  on  the  in- 


ventory  arc  be  V  units.  We  consider  two  alternatives  for  adjusting  the 
flow:  (1)  to  adjust  the  flow  such  that  the  flow  on  the  production  arc  van¬ 
ishes  and  that  on  the  inventory  arc  becomes  U+V,  or  (2)  to  adjust  the  flow 
such  that  the  flow  on  the  inventory  arc  vanishes  and  that  on  the  production 
arc  becomes  U+V. 

The  flow  adjustment  is  accomplished  by  creating  an  adjusting-flow  pat¬ 
tern  which  is  itself  feasible  at  every  node.  The  cost  associated  with  the 
pattern  will  be  concave,  thus  one  of  the  two  alternatives  will  yield  a  cost 
less  than  or  equal  to  the  original  cost. 

The  operation  of  adjusting  the  flow  at  a  given  node  which  violates  the 
WW  property  is  called  "operation-WWA" .  In  addition  to  yielding  a  WW-flow  at 
that  node,  the  operation  must  ensure  that  the  adjusted  flow  is  feasible,  and 
that  the  nodes  which  satisfied  the  WW  property  before  the  operation  must  con¬ 
tinue  to  satisfy  the  property  after  the  operation. 

Suppose  we  are  adjusting  the  flow  at  distributor  node  J.  The  adjusting 
**-  flow,  since-  it-  must  be  feasible,  should  correspond  to  a  "generalized  cycle- 
of  flow".  It  will  be  necessary  to  identify  not  only  the  set  of  arcs  in  this 
generalized  cycle,  but  also  the  direction  of  adjusting  flow  in  each  arc.  We 
construct,  to  serve  this  purpose,  two  subgraphs,  T^(J)  and  r2(J) ,  as  a  part 
of  operation  WWA.  r^(J)  contains  all  the  arcs  in  the  network  whose  flows 
contribute  to  the  flow  on  the  production  arc  at  J  and  similarly,  I*2(J)  con¬ 
tains  the  arcs  whose  flows  contribute  to  the  flow  on  the  inventory  arc  at  J. 
The  two  subgraphs  need  not  be  disjoint. 

In  a  given  feasible  flow  if  we  want  to  introduce  changes  in  the  inflow 
at  J,  corresponding  changes  have  to  be  introduced  on  arcs  in  T^(J)  and  r2(J) 
in  order  that  feasibility  of  flow  is  maintained.  To  satisfy  feasibility. 


the  changes  in  the  production  arc  and  the  entering  inventory  arc  at  J  are 
always  equal  in  magnitude  and  opposite  in  sign.  Changes  in  these  two  arc 
flows  induce  changes  in  T^(J)  and  r2(J)  respectively. 

The  example  network  of  Figure  2  will  be  used  to  illustrate  these 

notions.  ..Table  1  presents  a  set  of  arc  costs,  and  Figure  3  contains  a  fea- 

* 

sible  flow  which  violates  the  TO  property  at  node  23.  The  subgraphs,  T^(23) 
and  T2(23)  are  shown  in  Figure  4.  It  should  be  evident  from  this  figure 
that  the  changes  to  arc  flows  can  be  computed  "backward"  in  T^(23)  and  ^(23) 
starting  with  the  inventory  or  production  arc  and  using  the  arc  multipliers. 
Figure  5  shows  the  adjusted  flow. 

The  algorithm  for  constructing  a  TO  flow  is  given  in  Appendix  A.  The 
algorithm  makes  use  of  a  set  of  labels,  2^(m,  n)  defined  for  arcs 
(m,  n)  e  .  A  flow  change  of  A  on  arc  (J-T,  T)  -induces  a  change  of  - 

A^On,  n)  •  A  for  (m,  n)  e  r^(J)  and  -£2(m,  n)  *  A  for  (m,  n)  e  I*2(J).  Note 
that,  as  long  as  the  TO  property  is  satisfied  at  all  distributor  nodes  with1-?!.  .  •:.> 
indices  smaller  than  J,  if  an  arc  (m,  n)  is  in  both  r^(J)  and  I*2(J)  then 
i^(ra,  n)  ■  £2(m,  n).  Thus,  the  net  flow  change  for  such  an  arc  will  be  zero. 
Also,  the  adjusted  flow  cost  is  computed  by:  ' 

F(A)  -  l  C,  *(X(m,  n)  +  (£,(m,  n)  -  *2(m,  n))A) 

(m,  n)  e  rx(j)  u  r2( J)  °3’  n) 

Thus,  F(0)  is  tie  current  cost  for  arcs  in  the  two  subgraphs. 

Proof  of  Theorem  IN 

Let  J  be  a  distributor  node  such  that  the  TO  property  is  violated  by  J, 


and  is  satisfied  by  every  distributor  node  less  than  J.  Then  operation  TOA, 
when  applied  at  J  adjusts  the  flow  such  that: 


TABLE  1 


m 


Legend 

o - o  r2<23) 

o - o  r1<23) 


Figure  4.  Subgraphs  for  Flow  of  Figure  3 


1.  In  the  adjusted  flow,  the  WW  property  holds  at  node  J  and 
also  at  all  distributor  nodes  with  index  less  than  J. 

2.  The  adjusted  flow  at  nodes  greater  than  J  is  the  same  as  the 
flow  at  those  nodes  before  the  application  of  the  operation. 

3.  The  adjusted  flow  is  feasible. 

4.  The  cost  of  the  adjusted  flow  is  less  than  or  equal  to  the 
original  cost. 

The  first  two  of  the  above  four  properties  are  obvious  and  hence  we  omit 
a  formal  proof  for  those  properties. 

The  third  property  follows  from  the  fact  that  the  adjusting-flow  in  sub¬ 
graph  T^(J)  u  ^(J)  is  feasible.  Therefore,  the  adjusted  flow  which  is  the 
sum  of  the  original  and  the  adjusting-flows  should  be  feasible  also. 

To  establish  the  fourth  property,  we  compute  the  cost  changes  in 
r^(J)  u  I*2 (J)  and  show  that  the  total  cost  along  those  arcs  either  decreases 
or  remains  the  same  as  the  original  cost. 

First,  note  that  C,  % (X)  is  a  concave  function  of  X,  therefore 

tm,  n; 

C,  . (k  +  y)  is  a  concave  function  of  y  for  fixed  k.  Thus  F(A)  is  a  sum 
lm,  n; 

of  concave  functions  and  therefore  also  concave.  From  the  properties  of 
concave  function,  we  have: 

F(l)  >  F(0)  =>  F(-U)  <  F(0)  and 

F(l)  <  F(0)  =>  F(V)  S  F(0) . 

It  follows  that  the  cost  of  the  adjusted  flow  will  not  exceed  the  cost  of 
the  original  flow.  Q.E.D. 


THE  NESTED  SCHEDULE  PROPERTY 


A  nested  schedule  property  has  been  studied  by  Love  [2]  for  pure  series 
systems  and  by  Crowston  and  Wagner  [1]  for  pure  assembly  and  general  acyclic 
systems.  For  acyclic  systems,  a  nested  schedule  is  a  production  schedule  for 
which  production  at  any  stage  i  i  L  and  any  period  t  is  always  accompanied 
by  production  in  period  t  at  some  immediate  successor  of  stage  i. 

Crowston  and  Wagner  limited  their  discussion  to  the  case  of  unit  inter¬ 
stage  multipliers  and  linear  inventory  holding  costs.  We  allow  general  in¬ 
teger  valued  multipliers  and  concave  inventory  costs,  and  offer  a  network- 
based  proof. 


Theorem  2.  An  N-stage,  T-period  problem,  P,  has  at  least  one  optimal  solution 
that  simultaneously  satisfies  the  Wagner-Whitin  and  nested  sched¬ 
ules  properties  provided: 

.  .  1.  .‘.The  production  costs  for.  every,  stage-period  are  non-negative, 

concave,  and  zero  for  zero  production. 

2*  Within  a  stage  the  production  costs  are  nonincreasing  in 
time.  That  is 


X(X)  =  2(X)  * 


T(X) 


for  all  i 

3.  The  inventory  holding  cost  for  a  stage-period  is  non-negative 
and  concave  with  respect  to  the  end-of-the-period  inventory, 
and  zero  when  the  end-of-the-period  inventory  is  zero. 

4.  The  Inventory  holding  costs  are  non-decreasing  in  stages, 
and,  additionally,  satisfy  the  following  property: 


a 


i 

•  m 

£ 

1 


C 


fci 


=  I  H 

j  e  P(i) 


t  (Kji  X> 


for  i  i  F,  all  t. 

This  theorem  can  also  be  restated  in  terms  of  the  network  model: 


Theorem  2N:  Under  conditions  (l)-(4)  of  Theorem  2,  there  is  at  least  one 
optimal  flow  in  the  network  which  satisfies  the  Wager-Whitin 
property  and  in  addition,  whenever  there  is  positive  flow  in 
a  production  arc  (J-T,  J) ,  there  is  positive  flow  on  at  least 
one  tranfer  arc  originating  at  J. 

Our  proof  of  this  theorem  will  be  similar  to  the  proof  of  Theorem  IN;  a  flow 
adjustment  procedure  will  be  described  which  yields  the  desired  characteris¬ 
tic  without  increasing  total  cost. 

Initially  we  adjust  the  flow,  if  needed,  to  satisfy  the  WW  property  by 
applying  Algorithm  1.  At  this  phase  in  our  procedure,  if  there  are  some 
distributor  nodes  that  do  not  satisfy  the  nesting  property,  then  to  enforce 
nesting  we  apply  a  flow  adjustment  operation,  operation  USA.  However,  while 
enforcing  one  desirable  property  ,  nesting,  we  should  not  disturb  the  other 
desirable  property,  namely  the  WW  property,  which  the  flow  already  satisfies. 
Therefore,  in  case  the  WW  property  is  distrubed  by  operation  NSA,  it  will  be 
restored  by  another  operation,  operation  NSB. 

We  consider  nodes  from  2NT  to  0  consecutively  in  decreasing  order  of  in¬ 
dex;  when  a  node  that  does  not  satisfy  nesting  is  encountered  we  apply  first 
NSA  at  that  node  to  satisfy  nesting,  and  then  NSB,  if  needed  to  recover  the 
WW  property.  When  we  finish  considering  node  0,  we  will  have  a  nested,  WW 
flow  at  hand. 

Suppose  we  find  in  applying  the  above  procedure  that  we  need  to  apply 


L‘ 
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operation  MSA  at  node  J.  Since  J  violates  nesting,  J  must  have  incoming  flow 
on  its  production  arc  and  must  not  have  outgoing  flow  on  any  of  its  transfer 
arcs.  Also,  since  we  are  considering  nodes  in  the  descending  order,  all 
nodes  greater  than  J  must  satisfy  nesting,  and  the  WW  property.  Since  J 
violates  nesting,  and  since  the  flow  is  feasible,  J  cannot  belong  to  the 
first  period  nor  to  the  last  period  of  the  planning  horizon,  but  to  some  in¬ 
termediate  period.  Therefore,  there  exists  a  future  period  t  such  that  pro¬ 
duction  can  be  shifted  from  node  J  to  node  Q,  where  Q  is  the  distributor 
node  for  the  stage  under  consideration  in  period  t. 

Such  a  shift  of  production-flow  from  J  to  Q  necessitates  flow  changes 
in  several  ares*  The  arcs,  along  which  changes  are  to  be  made  are  contained 
in  subgraph  o(J) ,  the  construction  of  which  is  described  in  Algorithm  2.1.  ' 
The  proposed  shift  in  the  production  from  J  to  Q  is  obtained  by  sending  an 
adjusting,  flow  around  the  subgraph  o(J)-  The  adjusting  flow  is  so  computed 
that  the  resulting  flow  satisfies  the  following  properties:  (i)  the  node  J 
as  well  as  all  distributor  nodes  greater  than  J  satisfy  nesting  and  WW  prop¬ 
erties,  (ii)  all  nodes  satisfy  the  feasibility  conditions,  (iii)  the  cost 

/ 

of..  the_  resultant.,  flow.  Is  .not  ..greater  than  that  of  the  original  .flow.  . 

However,  in  the  resulting  flow,  some  nodes  with  indices  less  than  J 
may  not  satisfy  the  WW  property  after  changing  the  flow.  Therefore,  ve 
consider  nodes  from  1  through  J-l-  consecutively,  applying  operation  WWA  at 
those  nodes  where  WW  property  is  di9trubed  by  operation  NSA.  Such  consecu¬ 
tive  application  of  operation  WWA  from  1  through  J-l  is  referred  to  as  ap¬ 
plication  of  operation  MSB  at  J.  As  the  operation  NSB  does  not  distrub  the 
flow  beyond  node  J-l,  after  its  application  at  J,  the  WW  property  holds 
throughout  the  network.  The  complete  statement  of  the  algorithm  for  ob¬ 
taining  a  nested  schedule  is  given  in  Appendix  B. 


To  illustrate  the  procedure  for  obtaining  a  nested  schedule,  consider 
again  the  feasible  flow  shown  in  Figure  5.  This  flow  satisfies  the  WW  prop¬ 
erty,  but  the  nested  schedules  property  is  violated  at  node  23.  Operation 
NSA  Identifies  the  next  period  in  which  there  is  a  positive  flow  on  a  trans¬ 
fer  arc  out  of  this  stage  (3),  and  then  shifts  production  to  that  period. 
Shifting  the  production  means  that  flows  on  the  transfer  arcs  must  also  be 
shifted,  resulting  in  (possibly  additional)  inventory  carrying  at  the  pred¬ 
ecessor  stages.  Figure  6  shows  the  subgraph  o(23)  associated  with  the  nec¬ 
essary  flow  changes. 

After  this  flow  change,  the  WW  property  will  be  violated  at  node  16. 
Applying  operation  WWA  at  node  16,  we  find  that  node  7  violates  nesting. 
After  applying  operation  NSA  at  node  7,.  the  desired  result  is  achieved. 

This  flow  is  shown  in  Figure  7. 

Proof  of  Theorem  2N 

.  .  .In.  a.. given  feasible  flow: satisfying,  the  WW  property  let  J  be  a  distri¬ 
butor  node  that  does  not  satisfy  the  nesting  property.  Then  application  of 
operation  NSA  at  J  results  in  an  adjusted-flow  which  (1)  is  feasible,  (2) 
satisfies  WW  property  at  distributor  nodes  greater  than  or  equal  to  J, 

(3)  satisfies  nesting  property  at  J  and  at  all  distributor  nodes  greater 
than  J,  (4)  has  a  cost  less  than  or  equal  to  the  cost  of  flow  prior  to  ad¬ 
justment. 

Of  these  four  properties,  the  first  three  are  easy  to  see;  the  last 
property  is  not  so  obvious  and  hence  we  will  offer  here  a  proof  for  it. 

Note  that  in  the  operation  NSA  we  shift  production  in  the  stage  under 
consideration  to  a  future  period.  As  the  production  costs  are  non-increasln 
in  time,  the  production  costs  are  not  going  to  increase  as  a  result  of  this 
adjustment. 


The  inventory  costs  remain  to  be  considered,  however.  In  the  subgraph 
o(J) ,  for  every  inventory  arc  of  stage  i(J)  there  is  one  corresponding  in¬ 
ventory  arc  in  the  same  period  for  each  immediate  predecessor  of  i(J).  The 
flow  in  the  inventory  arcs  of  i(J)  decreases  by  U  units  and  that  in  the  in¬ 
ventory  arcs  of  the  predecessors  increases  by  K(.)U  where  K(.)  refers  to 
the  interstage  multiplier  between  i(J)  and  the  predecessor  stage  under  con¬ 
sideration. 

Consider  a  particular  inventory  arc  (J+K,  J+K+l).  This  arc  represents 
inventory  at  stage  j  =  i(J)  and  period  t  »  t(J+K).  Because  the  inventory 
costs  are  concave,  nondecreasing,  we  have: 


V*wV  -  V0)  4  W  +  KgjV  -  W  v  8-  '  a"d 


j  t  S(g) 


Also,  condition  4  of  Theorem  2  gives 


VV  %  l  Hj>  v(w- 

Because  Hjt(0)  =  0  for  all  j  and  t,  it  follows  that 


W  -  V0)  *  g  l  F(j)  <V<KgjV  -  v(0)} 


l 

g  C  rci) 


KgjV 


VV1 


Thus,  the  inventory,  cost  reduction  in  each  period  at  stage. i(J)  exceeds  the 
inventory  cost  increase  at  the  predecessor  stages.  Hence,  operation  NSA 


does  not  Increase  total  costs.  Q.E.D. 


CONCLUSION 


We  have  given  constructive  proofs  of  the  Wagner-Whitin  and  Nested  Sched 
ule  properties  for  optimal  solutions  to  acyclic  production  systems  with  gen¬ 
eral  integer  multipliers  and  concave  costs  of  production  and  inventory. 

These  proofs  are  important  for  two  reasons.  First,  since  they  are  network- 
based,  they  provide  greater  intuition  and  insight  into  this  complex  problem. 
Second,  because  they  are  constructive,  they  provide  a  simple,  efficient  meth 
od  for  obtaining  a  local  optimum,  given  any  trial  solution. 

These  results  can  be  extended  in  an  obvious  way  to  allow  backlogging  at 
the  final  stages.  The  existance  of  initial  inventories  is  difficult  to  for¬ 
malize,  but  only  affects  these  results  in  the  early  periods,  until  the  ini¬ 
tial  inventories  are  consummed.  The  case  of  constant,  identical  lead  times 
is  trivial;  proving  these  properties  with  nonidentical  lead  times  will  pro¬ 
bably  require  stricter  assumptions  on  the  cost  structures. 
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APPENDIX  A 


Algorithm  1:  (To  establish  the  Wagner-Whitin  Property) 
BEGIN 

FOR  J  =  1  to  2NT  DO 

IF  J  is  a  distributor  THEN 

IF  J  violates  WW  property  THEN 

•  Apply  operation-WWA  at  J  (See  Algorithm  1.1) 
END  IF 
ENDIF 
ENDFOR 

END 


Algorithm  1.1  (Operation-V>VA  at  a  given  node  J) 

BEGIN 

•  Construct  the  subgraph  T^(J)  (Algorithm  1.1.1) 

•  Compute  the  labels  L^(m)  for  each  node  of  T^(J) 
and  labels  l^(m,  n)  for  each  arc  of  T^(J) 
(Algorithm  1.1.2) 

•  Construct  the  subgraph  ^(J)  (Algorithm  1.1.3) 

•  Compute  labels  I^m)  for  each  node  of  ^(J)  and 
labels  ^(m,  n)  for  each  arc  of  r^CJ) 

(Algorithm  1.1. A). 

IF  F(l)  <  F(0)  THEN 

•  Adjust  the  flow  in  T^(J)  u  ^(J)  as  below: 
X(J-1,  J)  <■  0 

X(J-T,  J)  ♦  U+V 

X(m,  n)  •«-  X(m,  n)  +  (l^(m,  n)  -  n))V 

for  all  other  (m,  n)  e  ".(J)  u  r?(J) 


ELSE 


•  Adjust  tha  flow  in  T^(J)  u  r2(J)  as  below: 

X(J-1,  J)  «-  U+V 

X(J-T,  J)  0 

X(m,  n)  X(m,  n)  +  (l^m,  n)  -  l2Cn!,  n))(-U) 
for  all  other  (m,  n)  e  r^(J)  u  r2(J) 

ENDIF 

END 

Algorithm  1.1.1  (Construction  of  T^(J)) 

BEGIN 

•  Unraark  all  the  arcs  of  tha  network 

•  Mark  the  arc  (J-T,  J) 

•  FOR  P  -  J-T  to  0  DO 

•  IT  at  least  one  outgoing  arc  of  P  is  marked  AND 

there  is  at.' least  one  unmarked  incoming  arc 
with  a  positive  flow  THEN 

•  Mark  all  incoming  arcs  of  P  which  have  a 
positive  flow 

’  ENDIF 

•  endfor 

•  Include  all  marked  arcs  and  their  nodes  in  T^(J) 

END 
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Algorithm  1.1.2  (Construction  of  labels  in  T,(J)> 


BEGIN 

L^J)  «-  1 

FOR  b  *  J  to  1  DO 
IF  m  e  ^(J)  THEN 

IF  L^(m)  is  not  yet  computed  THEN 

FOR  each  outgoing  arc  of  m,  (m,  n) ,  in  r^(J)  DO 
IF  (m,  n)  is  a  transfer  arc  THEN 
l1(m,  n)  «*-  K(m,  n)  •  Lj(n) 

ELSE 

1  1 — ■  i 

^(b,  n)  «-  L-^n) 

ENDIF 

ENDFOR 

Ll(m)  "*  I  ^(m,  n) 

(m.njer^j) 

ENDIF 

ENDIF 

ENDFOR 

END 


Algorithm  1.1.3  (Construction  of  I\(J)) 

BEGIN 

•  Unmark  all  arcs  in  the  network 

•  Mark  the  inventory  arc  (J-l,  J) 

•  FOR  p  *  J-l  to  0  DO 

IF  there  is  at  least  one  outgoing  arc  of  P  which  is  marked 
AND  there  is  at  least  one  incoming  arc  of  P  with  a 
positive  flow  THEN 

*  Mark  all  incoming  arcs  of  P  which  have  a  positive  flow 
ENDIF 
ENDFOR 

Include  all  marked  arcs  in  T„(J). 


APPENDIX  B 


Algorithm  2:  To  establish  the  Nested  Schedule  Property 
BEGIN 

FOR  J  =  2NT  TO  1  DO 

IF  J  is  a  distributor  node  THEN 

IF  J  violates  nesting  property  THEN 

Apply  operation  NSA  at  node  J  (See  Algorithm  2.1) 
Apply  operation  NSB  at  node  J  (See  Algorithm  2.2) 
ENDIF 
ENDIF 
ENDFOR 

END 

Algorithm  2.1  (To  apply  operation  NSA  at  a  distributor  node  J) 
BEGIN 

•  Let  U  ^  X(J— T,  J) 

•"  Find  the  distributor  node  Q  in  stage  i(J)  such  that 

(i)  t(Q)  >  t(J) 

(ii)  At  least  one  outgoing  transfer  arc  of  Q  has 
a  positive  flow. 

(iii)  The  nodes  between  J  and  Q  do  not  have  flows 
on  their  outgoing  transfer  arcs. 

•  Include  inventory  arcs  (J,  J+l),  (J+l,  J+2),  ...  , 

;  (Q-l,  Q)  and  their 'end  nodes  in  the  subgraph  a(J). 

•  Let  the  adjusted  flow  in  arcs  (J,  J+l),  (J+l,  J+2), 

...  ,  (Q-l,  Q)  be  zero,  that  is 

X(J+n,  J+n+1)  =  0  for  n  *  0,  1,  ...  ,  Q-J-l 

•  Include  arc  (J-T,  J)  and  node  (J-T)  in  o(J) 


*.  Let  the  adjusted  flow  on  (J-T,  J)  be  zero,  that  is 
X(J-T,  J)  -  0 

•  Include  arc  (Q-T,  Q)  and  node  (Q-T)  in  a(J). 

>  Let  the  adjusted  flow  on  (Q-T,  Q)  be 

X(Q-T,  Q)  X(Q-T,  Q)  +  U 

-  IF  i(J)  i  F  THEN 

FOR  each  distributor  node  in  which  is  connected  to 
(J-T)  with  a  transfer  arc  (m,  J-T)  DO 

•  Find  the  distributor  node  r  in  stage  i(m)  such 
that  transfer  arc  (r,  Q-T)  exists. 

•  Include  transfer  arc  (a,  J-T)  in  o(J).  • 

•  Let  X(m,  J-T),  the  adjusted  flow  on  (m,  J-T), 
be  zero. 

-  •  Include  the  transfer  arc  (r,  Q-T)  in  o(J) 

•  Let  the  adjusted  flow  on  (r,  Q-T)  be 

X(r,  Q-T)  «•-  X(r ,  Q-T)  +  U  •  K(r,  Q-T) 

«- Include:  inventory  arcs  (a,  m+1) ,  (m+1,  m+2) ,  ... 
_.(r-l,  r)  and  the  nodes  a,  m+1,  ...  ,  r  in  o(J) 

•  •  Let- the.  adjusted  flow  along  the  inventory  arcs 

between  n  and  r  be 

X(m+n,  m+n+1)  X(a+n,  mfn+1)  +  U  •  K(m,  J-T) 
for  n  =  0,  1,  ...  ,  r-m-1 

ENDFOR 


Include  transfer  arcs  (0,  J-T)  and  (0,  Q-T)  and 
node  0  in  o(J) 

Let  X(0,  J-T)  *-  0 
Let  X(0,  Q-T)  -  X(0,  Q-T)  +  U 
Et?DIF 


Algorithm  2»3  (Application  of  operation  NSB  at  a  given  node  J) 
BEGIN 

FOR  r  »  0  to  J-l  DO 

IF  r  is  a  distributor  node  THEN 

IF  r  does  not  satisfy  the  WW  property  THEN 

Apply  operation  WWA  at  r  (see  Algorithm  1.1) 

ENDIF 

ENDIF 

ENDFOR 

END 
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